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We present a measurement of the Cabibbo-Kobayashi-Maskawa matrix element \Vcb\ based on 
a sample of about 53,700 -B" D*'^ l~vi decays observed by the BAS4R detector. We obtain the 
branching fraction averaged over £ = e, fj,, B{B'^ D'+rvi) = (4.90 ± 0.07(stat.)+Q 35(syst.))%. 
We measure the differential decay rate as a function of w, the relativistic boost 7 of the D*"*" in 
the rest frame. By extrapolating dF/dw to the kinematic limit ui — » 1, we extract the product 
of I Vet I and the axial form factor Ai{w = 1). We combine this measurement with a lattice QCD 
calculation of Aiiw = 1) to determine \Vcb\ = (38.7 ± 0.3(stat.) ± 1.7(syst.) + j-3(theory)) x 10"^. 

PACS numbers: 13.20.He, 12.15.Hh 



In the Standard Model of electroweak interactions, the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix describes 
the flavor mixing among quarks and determines the 
strength of CP violation by a single non-trivial weak 
phase. The CKM matrix element Vcb measures the 
weak coupling of the b to the c quark. In this Let- 
ter, we present measurements of the branching fraction 
B{B° D*+i-iyi) 1\ and \Vcb\- The rate for this weak 
decay is proportional to \Vcb\'^ and is influenced by strong 
interactions through form factors, which are not known a 
priori. In the limit of inflnite &-quark and c-quark masses, 
these form factors are determined by a single Isgur-Wise 
function 2]. The value of this function when the is 
at rest relative to the 5° has been computed for flnite c- 
and 6-quark masses using lattice QCD |3|. 

In this analysis, we measure the differential decay rate 
dT/dw, where w is the product of the four- velocities of 
the 5° and D*+, and corresponds to the relativistic boost 
7 of the D*~^ in the rest frame. We extrapolate the 
rate to the zero- recoil limit w=l, and use the theoretical 
result for the form factor there Q to extract \Vcb\- 

The analysis is based on a data sample of 79 fb~^ 
recorded on the T(4S') resonance and 9.6 fb~^ recorded 
40 MeV below it, with the BABAR detector |3| at the 
PEP-II asymmetric-energy e+e" collider. We use sam- 
ples of GEANT Monte Carlo (MC) simulated events that 
correspond to about three times the data sample size. 

The momenta of charged particles are measured by a 
tracking system consisting of a five-layer silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH), op- 
erating in a 1.5-T solenoidal magnetic field. Charged 
particles of different masses are distinguished by their 
energy loss in the tracking devices and by a ring-imaging 
Cherenkov detector. Electromagnetic showers from elec- 
trons and photons are measured in a CsI(Tl) calorimeter. 
Muons are identified in a set of resistive plate chambers 
inserted in the iron fiux-return yoke of the magnet. 

We select events that contain a D*^ candidate and 
an oppositely charged electron or muon with momentum 
1.2 < pi < 2.4 GeV/c. (Unless exphcitly stated other- 
wise, momenta are measured in the T(4S') rest frame, 
which does not coincide with the laboratory frame, due 
to the boost of the PEP-II beams.) In this momentum 
range, the electron (muon) efficiency is about 90% (60%) 
and the hadron misidentification rate is typically 0.2% 
(2.0%). We select D*^ candidates in the momentum 



range 0.5 < po' < 2.5 GeV/c in the channel _D*+ 
£'°7r+, with the 13° decaying to A'^7r+, K^tt^tt^tt^ , or 
A'^7r+7r°. The charged hadrons of the D'-' candidate are 
fit to a common vertex and the candidate is rejected if 
the fit probability is less than 0.1%. We require the in- 
variant mass of the hadrons to be within 17 MeV/c^ of 
the 13° mass for the decays to only charged particles, and 
34 MeV/c^ for R-tt+tt'^ decays. For D° K-tt+tt^, we 
accept only candidates from portions of the Dalitz plot 
where the square of the decay amplitude, as determined 
by Ref. 5], is at least 10% of the maximum it attains 
anywhere in the plot. For the pion from £)*"*" decay, nf , 
the momentum in the laboratory frame must be less than 
450 MeV/c, and the transverse momentum greater than 
50 MeV/c. Finally, the lepton, 7r+, and D° are fit to 
a common vertex with a beam-spot constraint, and the 
probability for this fit is required to exceed 1%. 

In semileptonic decays, the presence of an undetected 
neutrino complicates the separation of the signal from 
background. We compute a kinematic variable with con- 
siderable power to reject background by determining, for 
each B-decay candidate, the cosine of the angle between 
the momentum of the and of the D*^£~ pair, under 
the assumption that only a massless neutrino is missing: 



cos t^BO D' 



Ml 

BO 



Ml, 



This quantity constrains the direction of the to lie 
along a cone whose axis is the direction of the D*~^i~ 
pair, but with an undetermined azimuthal angle about 
the cone's axis. The value of w varies with this azimuthal 
angle; we take the average of the minimum and maximum 
values as our estimator w for w. This results in a reso- 
lution of 0.04 on w. We divide the sample into 10 bins 
in w from 1.0 to 1.5, with the last bin extending to the 
kinematic limit of 1.504. 

The selected events are divided into six subsamples, 
corresponding to the two leptons and the three 13° de- 
cay modes. In addition to signal events, each subsample 
contains backgrounds from six different sources: com- 
binatoric (events from BB and continuum in which at 
least one of the hadrons assigned to the D*'^ does not 
originate from D*~^ decay); continuum (D*'^£~ combi- 
nations from e~^e^ cc); fake leptons (combined with 
a true D*~^); uncorrelated background {£ and Z3*+ pro- 
duced in the decay of two different B mesons); events 
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FIG. 1: Yields of on-resonance data (points) and the results of 
the fit (line) to the AM distribution, with contributions from 
continuum, fake-lepton, and combinatoric-D*'*' backgrounds 
summed over all w bins. 



from B D*+7ri' i^i decays ; and correlated background 
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events due to the processes B° 
and B° D*+Xc, i^Y. We estimate correlated 

background (which amounts to less than 0.5% of the se- 
lected candidates) from Monte Carlo simulation based on 
measured branching fractions |6] , while we determine all 
the others from the data. Except for the combinatoric 
background, all other background sources exhibit a peak 
in the AAf ~ M£).+ — A/^jo distribution, where Af£)*+ and 
M/jo are the measured D*^ and candidate masses. 

We determine the composition of the subsamples in 
each w bin in two steps. First we estimate the amount of 
combinatoric, continuum, and fake-lepton background by 
fitting the AA/ distributions in the range 0.139 < AA/ < 
0.165 G<N/(? simultaneously to three sets of events: data 
recorded on resonance, data taken below the T{AS) (thus 
containing only continuum background), and data in 
which tracks that fail very loose lepton-selection crite- 
ria are taken as surrogates for fake leptons. The distri- 
butions are fit with the sum of two Gaussian functions 
with a common mean and different widths to describe 
D*^ D'^irf decays and empirical functions, based on 
the simulation, for the combinatoric background. The 
four parameters of the Gaussian functions are common, 
while the fraction of peaking events and the parameters 
describing the combinatoric background differ for the sig- 
nal, off-peak, and fake-lepton samples. 

Since the AAf resolution depends on whether or not 
the 7r+ track is reconstructed only in the SVT or in 
the SVT and DCH, the fits are performed separately 
for these two classes of events. We rescale the number 
of continuum and fake-lepton events in the mass range 
0.143 < AA/ < 0.148 GeV/c^, based on the relative 
on- and off-resonance luminosity and measured hadron 
misidentification probabilities. In the subsequent analy- 
sis we fix the fraction of combinatoric, fake-lepton, and 
continuum events in each w bin to the values so obtained. 
Figure 1 shows the AAf fit results for the on-resonance 
data. 
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FIG. 2: Yields of on-resonance data (points) and the results 
of the fit (histograms) to the cos Obo^d'i distribution, summed 
over all w bins. 



In a second step, we fit the cosOgo jjtg distributions 
in the range —10 < cos^^o j^*^ < 5 and determine the 
signal contribution and the normalization of the uncor- 
related and B D*'^'iTl~vg backgrounds. Neglecting 
resolution effects, signal events meet the obvious con- 
straint \cos9bo^d*i\ < 1) while B ^ D*^-Kl~vg events 
extend below —1, and uncorrelated background events 
are spread over the entire range considered. 

We perform the fit separately for each w bin, with 
the individual shapes for the signal and for each of 
the six background sources taken from MC simulation, 
specific for each of the six subsamples. Signal events 
are generated with the form-factor parameterization of 
Ref. 0, tuned to the results from CLEO % Radia- 
tive decays {W D*+£-iya) are modeled by PHOTOS 
9] and treated as signal. B D**lv decays involving 
orbitally excited charm mesons are generated according 
to the 1SGW2 model and decays with nonresonant 
charm states are generated following the prescription in 
Ref. To reduce the sensitivity to statistical fluctu- 

ations we require that the ratio of _B — s- D*'^tt£~ Pg and 
of uncorrelated background to the signal be the same 
for all three Z?" decay modes and for the electron and 
muon samples. Fit results are shown in Fig. [3 In total, 
there are 70,822 events in the range | cos fl^o.D'^ I < 1-2. 
The average fraction of these events that are signal is 
(75.9 ± 0.3)%, where the error is only statistical. 

To extract \Vcb\, we compare the signal yields to the 
expected differential decay rate 



dr 

dw 



MD.+ fg{w)J'{wf\Vcb\^ 



where 



g{w) = Vw^ - 1(10 + 1 2 1 + 4 — — 

\ w + 1 (1 ~ r)^ 

is a phase-space factor, r = M£)*+ /M^q. We parameter- 
ize the form factor J^{'w) with a Taylor expansion: 

T{w) « ~ pI{w c{w - 1)2), 
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where we neglect terms of order greater than two in {w — 
1). We fit the data to determine JF(l)|ycb|, P:^ and c. 

Dispersion relations inspired by QCD can be used to 
constrain the shape of the form factor and reduce the 
number of parameters to be determined 0, ^3 ■ There- 
fore we consider also the parameterization proposed in 
Ref. , which relates J-'{w) to the axial- vector form fac- 
tor Aiiw) according to the following expression: 



1 — 2wr + 



,w — 1 



X i + Ri{wy 



w + 1 



l + (l-i?2H) 



Ul ^ 1 

1-r 



0.05(w 

\2 



1) 



where Ri{w) w i?i(l) - 0.12(u; - 1) 
i?2M « i?2(l) + 0.11(ui - 1) - 0.06(w - 1)^ and we use 
the values i?i(l) = 1.18 ±0.32 and i?2(l) = 0.71 ± 0.21 
measured by CLEO Using dispersion relations we 
express the ratio Ai{'w)/Ai{l) as a function of a single 
unknown parameter p^_^ : 



Aijw) 
A(l) 



1 - 8pXz+ i53pX - 15)z2 - {231pX - 91)^^ 



where z = {y/w + 1- \/2)/(Vw + 1 + V2). It must be 
noted that, for w — s- 1, Ai{w) ^(w), so we expect 

A(l)«^(l). 

We perform a least-squares fit of the sum of the ob- 
served signal plus background yields to the expected yield 
in the ten bins in w . We define for each of the six data 
subsamples 



10 



bk 



i\2 



= 1 ^data 



'bk 



.7 = 1 j 



where N^^^.^ is the number of observed events in the i*^ 
bin; N^^ and cr^j. are the number of estimated back- 
ground events and its error. The backgrounds are fixed 
to the estimated rates. The expected signal yield is 
calculated at each step of the minimization from the 
reweighted sum of N^jq simulated events. Each weight 
is the product of four weights, = W-'' wj^'\ 
The factors W^, W-'' do not vary during the minimiza- 
tion, while the terms , W^^'^ depend on parameters 
which are determined by the fit, and vary at each step of 
the minimization. 

The first factor W'~ accounts for relative normaliza- 
tion of the data and MC samples, and is common to 
all subsamples. W'~ depends on the total number of BB 
events, Ng^ = (85.9 ±0.9) x 10*^, on the fraction of B'^B^ 
events, /oo = 0.489 ±0.012 "6^, on the branching fraction 
B{D*+ D°TT+) = 0.677 ± 0.005 and on the 5° 
lifetime — 1.536 ± 0.014 ps 0. Wj'* accounts for dif- 
ferences in reconstruction and particle-identification ef- 
ficiencies predicted by the Monte Carlo simulation and 



measured with data, as a function of particle momen- 
tum. Only the tt^ tracking efficiency varies significantly 
with w. 

The weight accounts for potential small differ- 
ences in efficiencies for the six data subsamples and al- 
lows for adjustments of the D° branching fractions, prop- 
erly dealing with the correlated systematic uncertain- 
ties. It is the product of several scale factors that are 
floating parameters in the fit, each constrained to an 
expected value with a corresponding experimental er- 
ror. For instance, to account for the uncertainty in the 
multiplicity-dependent tracking efficiency, we introduce a 
factor W^^^. = 1 + NtrkStrk, where Ntrk is the number of 
charged tracks in the D*^t~ candidates in each sample 
and 5trk is constrained to zero within the estimated un- 
certainty in the single-track efficiency: ±0.8%. Similarly, 
correction factors are introduced to adjust lepton, kaon, 
and 7r° efliciencies, and D'^ branching fractions, taking 
into account correlations. 

The fourth factor, wj^'^, adjusts the fitted decay dis- 
tribution relative to the one used in the generation of the 
MC events. This term depends on \ Vcb\ and on the shape 
parameters. It is a function of w and is determined for 
each simulated event at each step of the fit. 

Figure 13 (top) compares the observed signal and back- 
ground yields, summed over all six subsamples, with the 
result of the fit. Figure 01 (bottom) illustrates the extrap- 
olation to w = 1 for the two form-factor parameteriza- 
tions. The numerical values obtained for the two different 
form-factor parameterizations are listed in Tabled For 
both fits, the per degree of freedom is satisfactory, 
and the scale factors introduced to allow adjustments of 
the efficiencies and branching fractions deviate from their 
default values by less than one standard deviation. 



TABLE L Results of the fits to dT/dw for the two parameter- 
izations of the form factor. The errors stated include statis- 
tical error of the data and MC as well as uncertainties due to 
tracking, particle identification, and branching fractions 
that are directly assessed in the fit procedure. 



Ai{l)\Vcb\ X 10" 



xVndf 



T 35.0 ± 0.9 0.95 ±0.09 0.54 ±0.17 67/57 
Ai 35.5 ±0.8 1.29 ±0.03 - 69/58 



In Table II we present a summary of the statistical and 
systematic uncertainties. From the fit to the w distribu- 
tion we obtain errors that combine the statistical error 
with systematic errors introduced by the uncertainties in 
scale factors. We separate the various contributions in 
the following way: first, we extract the statistical errors 
by fixing all scale factors to their fitted values. The sys- 
tematic errors due to the uncertainties in a given scale 
factor is extracted from a separate fit in which this scale 
factor is fixed. We take the square root of reduction in 
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in 




1 1.1 1.2 1.3 1.4 1.5 

w 

FIG. 3: Results of the fit as a function of w compared to 
data. Top: the observed w distribution (points) compared to 
the fit result; signal and background contributions are indi- 
cated using the same shading as in Fig|5| Bottom: the form 
factor parameterizations with fitted parameters compared to 
the background- and efficiency-corrected data. The solid (dot- 
ted) line corresponds to the Ai{w) (T{w)) parametrization, 
and is to be compared to the filled (open) data points. 

the square of the fit errors as a measure of the contribu- 
tion of the particular scale factor to the overall error in 
the fit parameters. 

We then assess the individual contributions to the sys- 
tematic error due to other input quantities by varying 
their values by their estimated uncertainties and adding 
in quadrature the resulting changes to the fit parame- 
ters. The uncertainties in the lifetime tbo, the T{AS) 
and Z3*+ branching fractions, and overall normalization 
are independent of w and thus do not affect the shape 
of the form factor. The uncertainty introduced by the 
vertex reconstruction is common to all samples and in- 
dependent of w. It is determined by comparing the event 
samples with and without cuts on the vertex probability. 
The error induced by the cut on the decay amplitude for 
the K^tt^tt'^ decay is determined by varying that cut. 

A major source of uncertainty is the reconstruction ef- 
ficiency for the low-momentum pion from the D*~^ decay, 
since it is highly correlated with the D*"*" momentum and 
thereby with w. We determine the tracking efficiency for 
high-momentum tracks comparing the independent infor- 
mation from SVT and DCH. We compute the efficiency 
for low-momentum tracks reconstructed in the SVT alone 
from the angular distribution of the 7r+ in the D*'^ rest 
frame. We use a large set of D*+ -C'°7r+, D° K~Tr+ 
decays selected from generic hadronic events. For fixed 
values of the D*"*" momentum, we compare the observed 
angular distribution to the one expected for the decay of 
a vector meson to two pseudoscalar mesons. We define 



the relative efficiency as the ratio of the observed to the 
expected distribution and parameterize its dependence 
on the laboratory momentum of the 7r+. The study is 
performed in several bins of the polar angle of the detec- 
tor. We perform the measurement in the data and in the 
simulation, and we find that the functions parameterizing 
the efficiency are consistent within the statistical errors. 
To assess the systematic uncertainty on |Vc&|, we vary 
the parameters of the efficiency function by their uncer- 
tainty, including correlations. We add in quadrature the 
uncertainty in the absolute scale, as determined using 
high-momentum tracks reconstructed in both the SVT 
and the DCH. We obtain a systematic error of ±1.1% on 
\V.b\. 



TABLE II: Summary of uncertainties. 



Source of Uncertainty S{Ai 




l)(%) SpX 


5B{%) 


Data and MC statistics 


0.7 


0.03 


1.4 


B{D" ^ K--K+) 


1.1 




2.2 


B{D° ^ K-TT+TT-TT + ) 


0.4 




0.8 


B{D" K--r+t:°) 


0.5 




1.0 


Particles identification 


1.1 




2.2 


Tracking & tt" reconstr. 


1.3 




2.6 


Partial Sum 


2.2 


0.03 


4.5 


B" lifetime 


0.5 






Number of BB 


0.6 




1.2 


B{D*+ ^ D°Tv+) 


0.4 




0.7 


S(r(4S) B°B") 


1.2 




2.5 


D*'^£~ vertex efficiency 


0.5 




1.0 


TTs efficiency 


1.1 


0.01 


1.9 


D*nlv sample composition 


1.8 


0.06 


2.0 


B momentum 


0.3 




0.7 


Radiative corrections 


0.2 


0.01 


0.4 


cos 6 go I & i5 fit method 


0.8 


0.02 


1.6 


i?i(l) and i?2(l) 


+2.9 

-2.6 


0.26 


+3.9 
-3.3 


Total Error 


-1-4.6 
-4.4 


0.27 


-1-7.4 
-7.1 



The largest error in the background subtraction is due 
to the uncertainty in the composition and form factors 
of the D*~^Tr£~Pe decays. We consider twelve different 
D*^TT states, narrow and wide, as well as nonresonant 
D*^TT. To assess the impact of these decays on the fit 
we repeat the analysis assuming that only one mode at 
a time populates the whole sample, and then take as the 
systematic error half the difference between the maxi- 
mum and minimum fitted parameters. 

We assess the effect of the uncertainty in the average 
momentum, as determined from a sample of fully re- 
constructed hadronic B decays on the fit results. We take 
into account an uncertainty of ±30% in the emission rate 
of the radiative photons predicted by PHOTOS j^. 

We also assess the impact of changes in the bin size on 
the fits to the cos p. ^ and w distributions. 

There are several uncertainties related to the form fac- 
tors and their parameterization. The form factor ratios 
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i?i and i?2 affect the lepton momentum spectrum and 
thus the differential decay rate as a function of w, as well 
as the fraction of events satisfying the lepton momen- 
tum requirements. We assess these effects by varying Ri 
and i?2 within the measurement errors taking into 
account their correlation. As a consistency check, we 
compare the measured momentum spectra of the D*~^ 
and leptons with the spectra expected from the fit re- 
sults. We find very good agreement for the Z?*+, but the 
lepton spectrum favors a larger value for Ri , though one 
consistent with the available measurement. 

If we fit separately e and samples, we find exactly 
the same value for . The values of y^ilKhl , (35.8 ± 
0.5) X 10-3 and (35.0 ±0.5) x 10'^ respectively, differ by 
1.2 standard deviation. 

The value of c, given in Table shows that the 
data disfavor a purely linear dependence of on w, 
by almost three standard deviations. The fits for the 
two different parameterizations of the w dependence of 
the form factors are consistent a.t w = 1. We choose 
Ai{l)\Vcb\ = (35.5 ± 0.3 ± 1.6) x lO^^, and = 
1.29 ± 0.03 ± 0.27, where the errors listed refer to the 
statistical, and the systematic uncertainties. The corre- 
lation between ^i(l)|T46| and is 0.56, taking into 
account statistical and systematic errors. A recent lat- 
tice calculation 3] (including a QED correction of 0.7%) 
gives Ai{l) = = 0.919+{]:[]^^, with which we obtain 

\Vcb\ = (38.7 ± 0.3 ± 1.7 li j) X 10-^ 

where the first error is statistical, the second is system- 
atic, and the third reflects the uncertainty in ^i(l). Inte- 
grating over the fitted w distribution these parameters 
result in the branching fraction B{B'^ D*~^£~Di) — 
(4.90 ± Omt"o.lt)%, where the errors are the statistical 
and systematic uncertainties. 

In summary, we have measured the CKM parame- 
ter \Vcb\ and the exclusive branching fraction for 
D*~^£~Di with high precision. The result for \Vcb\ is con- 
sistent with another BABAR measurement based on lep- 
ton and hadron spectra from inclusive semileptonic B- 
meson decays [3, \Vcb\ = (41.4±0.4(stat.)±0.4(exp.)± 
0.6(theory)) x lO""^. The results for | Vch| and the branch- 
ing fraction are also consistent with earlier measure- 
ments |l4j based on the technique employed here, except 
for those from the CLEO experiment |l5|. 
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